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We systematically construct tetraquark currents of I'^ and classify them into types 

A (antisymmetric), S (symmetric) and M (mixed), based on flavor symmetries of diquarks and 
antidiquarks composing the tetra quark currents. We use tetraquark currents of type M to perform 
QCD sum rule analyses, and find a tetraquark current with quark contents qsqs {q = u or 
d) leading to a mass of 1.44 ± 0.08 GeV consistent with the ai(1420) state recently observed by 
the COMPASS collaboration. Our results support tetraquark explanations for both ai(1420) and 
/i(1420), assuming that they are isospin partners. We also study their possible decay patterns. As 
tetraquark candidates, the possible decay modes of ai(1420) are S-wsve ai(1420) ^ K*{892)K and 
P-wave ai(1420) —>■ /o(980)7r while the possible decay patterns of /i(1420) are P-wave /i(1420) —>■ 

K*{892)K and P-wave /i(1420) —>■ ao(980)7r. We speculate that ai(1420) is partly responsible for 
the large isospin violation in the 77(1405) —>■ /o(980)7ro decay mode which is reported by BESIII 
collaboration in the J/t/) —>■ 737r process. 

PACS numbers: 12.39.Mk, 11.40.-q, 12.38.Lg 
Keywords: tetraquark, axial-vector meson, QCD sum rule 


I. INTRODUCTION 

Recently, the COMPASS collaboration at CERN observed a narrow = 1++ signal in the /o(980)7r channel, 
and identified a new ai state with mass MeV and width 153)')23 MeV [IHi]. Including ai(1260), ai(1420), 

ai(I640), ai(1930), ai(2095) and ai(2270), there are as many as six ai states of quantum numbers l‘^ 

(see Ref. @ and references therein), which is much richer compared with conventional quark model predictions. 
Moreover, this new ai(1420) state was observed in the /o(980)7r channel, suggesting that ai(1420) has a large ss 
component, since /o(980) is usually interpreted as a KK molecule, diquark-antidiquark tetraquark, or other models 
with an ss component (THlsjl. In particular. Refs. (isl - fTsjl used the diquark-antidiquark model to interpret the light 
scalar meson spectrum. Accordingly, ai(1420) can be an exotic multiquark state, which makes it quite an interesting 
subject. 

To date, there are only a few theoretical studies of ai(1420). In Ref. [T^, it was interpreted as an axial-vector 
two-quark-tetraquark mixed state using QCD sum rule methods, but the analysis is incomplete because it did not 
include the interference terms in the calculation of the mixed correlator. Ref. [13 interpreted the ai(I450) as a 
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dynamical effect due to the singularity in the triangle diagrams formed by the processes ai(1260) —>■ K*K, K* Kn 
and KK /o(980)(+c.c). It was also briefly discussed using lattice QCD in Ref. [T^ . 

In this paper we shall study the ai(1420) state in the framework of QCD sum rules, which has proven to be a 
successful and powerful nonperturbative method over the past few decades [il,[2^. We shall systematically construct 
local interpolating tetraquark currents of and classify them into types A, S and M, respectively 

based on antisymmetric, symmetric, and mixed flavor symmetries of diquarks and antidiquarks. Tetraquark currents 
of types A and S have been investigated in Ref. (2lj |. while in this paper we shall use tetraquark currents of type M 
to perform QCD sum rule analyses. We shall find a tetraquark current with quark contents qsqs {q_= u oi d), which 
leads to a mass result consistent with the ai(1420) state observed by the COMPASS collaboration 0. Possible decay 
patterns based on this current will also be studied. 

This paper is organized as follows. In Sec. im we systematically construct tetraquark currents of I'^ = 1 1++. 

Then in Sec. IIIII we use these currents of type M to perform QCD sum rule analyses. In Sec. IIVI we summarize our 
results and discuss possible ai(1420) decay patterns. 


II. INTERPOLATING FIELDS 


The flavor structure of light tetraquarks is 

3(g)3(g)3(g)3 = (3 ©6) 0 (3 © 6) (1) 

= (3^3) © (3 0 6) © (6 0 3) © (6 0 6) 

= (l © 8) © (8 © To) © (8 © 10 ) © (l © 8 © 27) , 

where the subscripts A, M 1 /M 2 and S denote that the diquarks and antidiquarks inside have antisymmetric, mixed- 
symmetric and symmetric flavor structures, respectively. 

The tetraquark currents of = I’*' have been constructed in Ref. [2l|. Now we need to take the charge-conjugation 
parity into account in order to construct tetraquark currents of The charge-conjugation transforma¬ 

tion changes diquarks into antidiquarks, and vice versa, while keeping their flavor symmetries unchanged. Therefore, 
tetraquark currents themselves can have definite charge-conjugation parities when the diquark and antidiquark fields 
inside have a symmetric flavor structure 6f (gg) 0 6f{qq) (S) or an antisymmetric flavor structure 3[{qq) 0 3f (gg) (A). 
These currents have been constructed in Ref. [2l|, i.e., there are two independent tetraquark currents of quantum 
numbers and type S: 

= qA'^<iB{dvii^i5<Cq’S' + ^cit^i5<Cq^) + q‘fCjf,j5q’h{qv^qW + dc^d‘D), ( 2 ) 

^/’s,2a. = q'f<CYqBiqc^t^''^5<Cq’0' - gcO-^xi-dsCgj)^) + q‘f “ qci''^qD') > 

and two independent tetraquark currents of quantum numbers = 1“''^ and type A: 

= qf + qfCj^lMq^cCq^D ” , (3) 

V’A.2a. = qf'<C'y''q%{q^afj,^-f5Cq'0' + + qf‘Ca^.^-f^q^iq^Y^Cq’S' + ■ 


In these expressions, q\(x) = [Ma(a;) ,da{x) ,So(a^)] denotes the flavor triplet quark field; A---D are flavor indices; 
a and b are color indices; C is the charge-conjugation operator; and the superscript T denotes the transpose of the 
Dirac indices. 

Tetraquark currents constructed using combinations of 3f 06f (Mi) and 6f 03f (M 2 ) can also have definite charge- 
conjugation parities (see Ref. for detailed discussions). Tetraquark currents of = I’*' and types M 1 /M 2 have 
been constructed in Ref. [Hj |: 

V'MiPm = {qA’<Cl5q’h){qllt,Cq’’J - q’Xj^Cq"^) , 

V'Mi,2m = (9A^C7^9s)(9n75Cg^^ + 475C^^), 

V'Mi.Sm = {qf'^lulbqB){q\<x^^<Cq'^i^ - q\a^^Cq^i^) , 

V'MiAm = (g(4^CCT^i.gB)(9A7i^75Cg^ + g^7j.75C^^), (4) 

= (9n^C7^9s)(9n75Cg^^-g^75Cg27’), 

V'M2,2m = (9A^C75gB)(gA7A*Cg^^+ 47^Cg^^), 

V'M2,3m = (9(4^Ccr^i^g^)(g(^7i,75Cg^^ - g^7j.75Cg5,^), 

V'M2,4m = {qf'Cj^^5q%){q^af,,,Cq’g' + . 
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We can use these currents to construct positive charge-conjugation parity currents = I'*'"'") 

V'M.i/i = = li ■ ■ • ; 4 , (5) 

as well as negative charge-conjugation parity currents {J^'^ = 1“' ) 

5 ^ 4 ; ' * ' ; 4 , ( 6 ) 

where we have denoted them as types M and M'. 

We now consider the isospin degree of freedom (see Fig. 1 of Ref. [22| and related discussions). There are two 
isospin triplets of type S, whose quark contents are 

OTOT(S) , qsqs{S) - 6f (g) 6f (S) , (7) 

one isospin triplet of type A, whose quark contents are 

qsqs{A) ^ 3f 0 3f (A), (8) 

and two isospin triplets of type M, whose quark contents are 

qqqq{M) , qsqs{M) ~ (3f 0 6f) 0 (6f 0 3f) (M) . (9) 

In these expressions q represents an up or down quark, and s represents a strange quark. 

With all these analyses, we can construct tetraquark currents of quantum numbers = 1“1++ and types 

S/A/M and collect them as follows. 

1. For the two isospin triplets belonging to 6/ 0 6y (S), there are altogether four independent tetraquark currents. 
Among them, two contain only light flavors, and the other two contain one ss quark pair: 

= i’s.ifj.iqqqq) ^ ulCdbiualfi'JbCd^ + Ub-ffj.-f5<Cd^)+ulCjf_,'y5db{uaCd[+UbCd^), 

V 2 bL = i’s. 2 fj.iqqqq) ^ulC'y''db{ua(Tf^^j5Cd[-UbCTf^^-fbCd^)+u^Caf_,^jbdbiual''<Cd^-Ubj''Cd^), (10) 

Vifj. = 'ips,ltj.{qsqs) U^CSf,(Ua7/x75Csf + UbJf,-f5Csl) + uf C7^75Sb(UaCsf + UbCsl) , 

vifi = i’s. 2 fj.iqsqs) ulC'j''Sb{Ua<7- UbCFfi^^b'Csl) + ulCcFfj,y^b,Sb{Ual''^sl - UbY’^s'^) . 

2. For the isospin triplet belonging to 3/0 3/ (A), there are two independent currents. They both contain one ss 
quark pair: 

Vlf, = IpAStiiqSqs) ulCsbiUa'Yfj.'ybCsl - UbJf,-f5Csl) + ulC'Yfj^'JbSbiUaCsl - UbCsl), (11) 

'n2^l = i’A,2tiiqsqs) ulCYsb{Uaaf_,^JbCsl + Ubafj,^j5Csl) +U^Caf_t^j5Sb{UaY‘Csl + UbYCs'^). 

Note that the two corresponding currents with quark contents udud, ipA.ifj.iqqqq) and '4’A,2fj.i,qqqq), both have 
isospin zero, as shown in Fig. 1 of Ref. 

3. For the two isospin triplets belonging to (3/06/) 0(6/03/) (M), there are eight independent currents. Among 
them, four contain only light flavors, and the other four contain one ss quark pair: 

r]f^l = 4’Ai,iti{qqqq) ^ u^Cj5db{ua^fj.Cd[ - Wby^Cd^) + u^C 7 ^dt,(ua 75 Cd^ - Ub-f 5 <Cd ^), 

= i’m, 2 ti{qqqq)ulCjf^dbiualb’Cd^ + UbjbCd^) + u'^Cj5db{ua^fj.Cd[ +Ubjf_,Cd^), 

= 'ipm,3tiiqqqq) u^Cj^-fbdb{uaCrfj,^Cd[ - UbU^^Cd^) + v^CcFfj,^db{ua'yulb^dl - Mb7i'75CdJ), 

= ’ipAiAt^iqqqq) u^Caf^^dbiuaj^j5<Cd[ + Ubj^-fbCd^) + v^C^ulbdbiuaCftiv'Cd^ + Ubcr^^CdJ), 

= V'M.i/i(gsp) ~ ulCjbSbiua'Yfj.Csl - Ub'jfj.Csl) + ulC'yf.Sbiua'lbCsl - Ub'YbCsl ), (12) 

r]^^l = 4 ’Ai, 2 ii{qsqs) ulCjf_,Sb{Ual5‘Csl +i2b75Cs^) + ulC'^bSbiua^fj.Csf + ltb 7 ^Cs^) , 

= i’m,3(i{qsqs) U^C-fulbSbiUaCTfiyCsl - ItbCT^yCs^) + U^<Caf_,ySb{Ua'yvl3^sl - UbJi^JbCsl) , 
r]^^l = 'ipm,itiiqsqs) ~ V^'Ca^uSbiUalvlb'Csl + Ub'^ulb'Cs'l) + U^C'^^-fbSbiUaCTt^^Csl + Ubcr^^Cs^). 
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In these expressions the quark content is not exactly correct, so we use instead of As an example, the 

current 77 ^^ contains quark content usus and so does not have 1=1. To be isovector, it should have quark content 
{usus — dsds), while its isoscalar partner should have quark content (usus + dsds). However, we do not study effects of 
isospin breaking in this paper, i.e., we work in the limit of SU(2) isospin symmetry and ignore the difference between 
up and down quarks, such as their masses and the quark condensates {uu) and {dd). Similarly, isospin-violating 
effects from instantons, which are important in the scalar channels [ 2 ^, are suppressed for the vector channel under 
consideration. Accordingly, the QCD sum rule results for these two currents with quark contents (usus — dsds) 
and {usus + dsds) are both the same as the result for 77 ^^ with quark contents usus and similarly for the other 
currents listed above. This suggests that we would obtain the same sum rule for an isovector tetraquark current of 
iGjPC _ 1 -X+- 1 - anft isoscalar partner of = 0 + 1 ++, which would consequently result in the same mass 

result for the relevant isovector state and its isoscalar partner. 

The tetraquark currents of types A/S and quantum numbers = 1“1++, rif^{qqqq), ilifj.{Qsqs), 

vifj.idsqs), r]f^{qsqs) and 7?^((7'S9s), have been used to perform QCD sum rule analyses in Ref. [2l|, and respectively 
result in similar masses, 1.51 - 1.57 GeV, 1.52 - 1.57 GeV, 1.56 - 1.62 GeV, 1.56 - 1.62 GeV, 1.57 - 1.63 GeV and 
1.57 — 1.63 GeV(see Sec. 5.2 and Fig. 4 of Ref. [2l| for detailed discussions, where all mass curves have a minimum 
around 1.5-1 .6 GeV against the threshold value sq)- One conclusion of Ref. [2l| is that these tetraquark currents 
couple to the ai(1640) state. Recently, a new ai(1420) state was observed, with mass 1414/^ijg MeV and width 153/^23 
MeV The masses of these two states are not far from each other, so that if a current couples to both of them and 
we still use a one-pole parametrization (see Eq. (11611 below and related discussion), a prediction between these two 
masses would be obtained for this single pole model. This may be the reason why the mass prediction 1.5-1 .6 GeV 
is obtained in Ref. [2l|. 

To better understand the properties of ai(1420), we need to differentiate it from ai(1640). To do this one can 
either adopt a two-pole parametrization, or use a current mainly coupling to ai(1420). The former is impractical 
because one needs detailed phenomenological models for such closely-spaced resonances, so in this paper we shall try 
the latter approach. Gonsidering that only tetraquark currents of types A and S were investigated in Ref. [2l|, we 
shall use tetraquark currents of type M, 77 )^ (i = 1 • • • 8), to perform QGD sum rule analyses and check whether such 
a current exists or not. We assume that they couple to the oi state of = 1“1++ through 

(0K|ai) = /M.zeM,* = l--- 8 > (13) 

where /m,z is the decay constant. 


III. QCD SUM RULE ANALYSIS 


We consider the following two-point correlation function 


= I j dW«"(o|r4(x)ji(o)|o) (14) 

= n(77^)(g,.-^) + n'(77^)^, 

q q 

in which (x) is an interpolating current carrying the same quantum numbers as the hadron state we want to study. 
Because Jfi{x) is not a conserved current, there are two different Lorentz structures in H^i/, n(( 7 ^) and n'(g^) related 
to spin -1 and spin -0 states, respectively. 

The two-point function Il^^{q'^) can be calculated in the QGD operator product expansion (OPE) up to certain 
order in the expansion, which is then matched with a hadronic parametrization to extract information about hadron 
properties. To do this, we express Eq. (HI at the hadron level as 


n(g^) 



Imn(s) 
s — q^ — ie 


(15) 


where we have used the form of the dispersion relation with a spectral function with s< denoting the physical threshold. 
We can write the imaginary part of Eq. m as 


Imn(s) = Hs - M{^){0\rj\n){n\r]'<\0). 

n 


(16) 


As usual, we adopt a parametrization of one-pole dominance for the ground state and a continuum contribution, but 
note that the masses of ai(1420) and ai(1640) are not far from each other so that it may be more reasonable to adopt 
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a two-pole parametrization, which is, however, impractical because one needs detailed phenomenological models for 
such closely-spaced resonances. After performing Borel transform at both the hadron and QCD levels, the two-point 
correlation function can be expressed as 

= - / e-"/^sinin(s)ds. (17) 

Finally, we assume that the contribution from the continuum states can be approximated well by the OPE spectral 
density above a threshold value sq (duality), and arrive at the sum rule relation which can be used to perform 
numerical analyses. Here we again use the current r]^ = tjjfj^{qsqs) as an example, whose quark contents are qsqs. 
We assume it couples to the ai state through Eq. (Il3l) . and the obtained sum rule relation is listed in Eq. (IA5I) . 
The results for other currents are shown in Appendix. We note that the Mathematica FEYNCALC package 
is used to calculate these OPEs. In these equations, there are dimension D = 3 quark condensates (qq) and (ss), 
D = A gluon condensate {g^GG), and D = 5 mixed condensates {gqaGq) and (gsaGs). The vacuum saturation for 
higher dimensional condensates are assumed as usual, such as {0\qqqq\0) ~ (0|g(7|0)(0|g(7|0). We have neglected the 
chirally suppressed contributions from current up and down quark masses because they are numerically insignificant. 
Moreover, we consider only leading-order contributions of as from the two-gluon condensate {{g^GG)) because the 
terms containing quark-related condensates are found to be significantly larger than those containing gluon-related 
condensates. 

To study the convergence of Eq. (IA5D , we use the following values for various condensates [2^ - 13^ : 

{qq) = -(0.240 ± 0.010)^ GeV^ , 

{ss) = -(0.8 ± 0.1) X (0.240 GeV)^ , 

= (0.48 ± 0.14) GeV^ , (18) 

{gsqcrGq) = -M^ x (qq) , 

(gsSaGs) = -Mq x (ss), 

= 0.8 GeV^ , 

ms(l GeV) = 125 ± 20 MeV. 

Note that there is a minus sign implicitly included in the definition of the coupling constant gs in this work. We find 
that the D = 6 and D = 8 terms are dominant power corrections, while the D = 10 and D = 12 terms are much 
smaller. Actually, the D = 6 and D = 8 terms in Eq. (|A5I) are mainly contributed by the condensates {qq){ss) and 
{qq){gsSaGs)/{ss){gsq<jGq), respectively. Accordingly, our first criterion is to require that the D = 10 and D = 12 
terms be less than 10%: 


Gonvergence (CVG) = | 


nM,5(oO,M|) 


I < 10%, 


(19) 


where H^®^ “‘^®“^(so. Mg) is the sum of the D = 10 and D = 12 terms. We show this in the left panel of Fig.[T] which 
shows that the OPE convergence improves with the increase of M^. This criterion has a limitation on the Borel mass 
that Mg > 1.1 GeV^. We note that this criterion gives almost no limitations if we assume Al^) to 

only contain the D = 12 terms, which implies that the contribution of the D = 12 terms is numerically small. 

Our second criterion is to require that the pole contribution be larger than 10% (see discussions below for the 
limitation 20%): 


Pole contribution (PC) = > 10% . (20) 

nM.5(oo,M|) - 

We note that the pole contribution is usually quite small in the multi-quark sum rule analyses due to the large powers 
of s in the spectral function. We show the variation of the pole contribution with respect to the Borel mass Ms in 
the right panel of Fig. [TJ when sq is chosen to be 2.5 GeV^. It shows that the PC decreases with the increase of Mb- 
This criterion has a limitation on the Borel mass that Mg <1.5 GeV^. Finally we obtain the working region of Borel 
mass 1.1 GeV^ < Mg <1.5 GeV^ for the current = ii}^^{qsqs) with the continuum threshold sq = 2.5 GeV^. 

Our final expression for the mass of the ai state is obtained via: 

11^,5 (so I Af^) 


( 21 ) 
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Borel Mass" [GeV"] 



FIG. 1: In the left panel we show CVG, as defined in Eq. m, as a function of the Borel mass Mb ■ In the right panel we show 
the variation of PC, as defined in Eq. PH). as a function of the Borel mass Mb- The current = '>lj^{qsqs) is used here and 
the threshold value is chosen to be so = 2.5 GeV^. 




FIG. 2: The mass calculated using the current = tpifiiqsqs), is shown with respect to the threshold value sq (left panel) 
for Mb = 1.1 (dotted), 1.3 (solid) and 1.5 GeV^ (dashed), and with respect to the Borel mass Mb (right panel) for so = 2.3 
(dotted), 2.5 (solid), and 2.7 GeV^ (dashed). The working region is 1.1 GeV^ < M% < 1.5 GeV^. 


in which Sq is the continuum threshold. To choose a reasonable value of Sq: show the variation of with respect 

to the threshold value sq in the left panel of Fig. [21 in a large region 1.5 GeV^ < sq < 3.5 GeV^. We find that the 
dependence of the mass curves with respect to the Borel parameter is very weak when the continuum threshold 
So is chosen to be around 2.5 GeV^, which is thus a reasonable value of sp to give a reliable mass prediction. 

The variation of with respect to the Borel mass Mb is shown in the right panel of Fig. [21 in a large region 0.5 
GeV^ < Mq <2.5 GeV^. The mass curves increase quickly with from 0.5 GeV^ to 1.0 GeV^, but they are quite 
stable against Mg as it continues increasing from 1 GeV^. This suggests that the limitation of the second criterion, 
Eq. (I2II1) . can be slightly modified to be 20%, and then the mass obtained is almost the same, but with a much narrower 
working region. Finally, we choose 2.3 GeV^ < Sq < 2.7 GeV^ and use the Borel window 1.1 GeV^ < Mg < 1.5 GeV^ 
as our working region resulting in the following numerical results 

Mai = 1.44 ± 0.08 GeV, (22) 

5 = (1.9 ± 0.5) X 10-3 GeV® , (23) 

where the central values correspond to sq = 2.5 GeV^ and Mg = 1.3 GeV^. The errors come from the uncertainties 
of so: Mg and the various parameters in Eq. (IPl) . The coupling constant /m .5 defined in Eq. (IT2)) gives the strength 
of the overlap between the interpolating current r]^ and the ai(1420) state. 

The sum rule using the current = ip^{qsqs) is similar to the previous sum rule obtained using the current 
r]^ = 'tpfl^{qsqs). The results are shown in Appendix PI and Fig. [3] We again choose 2.3 GeV^ < so < 2.7 GeV^ and 
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FIG. 3: The mass calculated using the current is shown with respect to the threshold value so (left panel) 

for M% = 1.1 (dotted), 1.3 (solid) and 1.5 GeV^ (dashed), and with respect to the Borel mass Mb (right panel) for so = 2.3 
(dotted), 2.5 (solid), and 2.7 GeV^ (dashed). The working region is 1.1 GeV^ < M% < 1.5 GeV^. 


use the interval 1.1 GeV^ < Mg <1.5 GeV^ as our working region, and obtain the following numerical results: 

M„, = 1.50 ± 0.08 GeV, (24) 

/m.6 = (2.6 ± 0.7) X 10-3 GeV® , (25) 

where the central values correspond to so = 2.5 GeV^ and Mg = 1.3 GeV^. 

The sum rules using the currents and r]^ = ^’^(9999) to larger masses around 1.6 GeV, as 

shown in Appendix El and Fig.Sl We choose 2.8 GeV^ < so < 3-2 GeV^ (we note that Sq should be larger than M^J 
and use the interval 1.2 GeV^ < Mg <1.8 GeV^ as our working region. We obtain the following numerical results 
for ? 7 (^: 


Mia, = 1.61 ±0.06 GeV, (26) 

= (3.0 ± 0.8) X 10-3 GeV® , (27) 

and the following numerical results for 77 ^: 

Ma, = 1.64±0.08 GeV, (28) 

/m ,2 = (4.2 ± 1.3) X 10-3 GeV® , (29) 


where the central values correspond to sq = 3.0 GeV^ and Mg = 1.5 GeV^. One notes that the central-value masses 
obtained in Eqs. (EZl) and (1291) for the non-strange qqqq tetraquarks are heavier than those in Eqs. (051) and (051) for 
the strange-flavor qsqs tetraquarks. This counter-intuitive behavior also appears in the scalar meson sector, where 
the isovector ao(1450) is a bit heavier than the strange isospinor ArQ(1430) [ 2 ^. However, it is important to note 
that the dominant sources of theoretical uncertainty in the QCD input parameters Eq. dni) are uncorrelated in the 
strange and non-strange cases, and hence we cannot rule out a near degeneracy from our mass predictions. 

The sum rules using the currents 77 ^ = ipif^{qqqq), 774 ^ = V' 4 ^(<?ggg), Vy, = V' 3 ^('?sgs) and = ip^^{qsqs) do not 
have reasonable working regions to give reliable mass results. We show these results in Appendix The sum rules 
using the former two currents lead to mass results roughly around 1.6 GeV, suggesting that they may couple to the 
ai(1640) state. The sum rules using the latter two currents lead to mass results around 1.8 GeV. 


IV. SUMMARY AND DISCUSSIONS 

In summary, we have systematically constructed tetraquark currents of I'^ = 1-1“''+. These currents can be 

classified into types A (anti-symmetric), S (symmetric) and M (mixed structure), based on flavor symmetries of 
diquarks and antidiquarks. Tetraquark currents of types A and S had been studied in Ref. [2lj, and in this paper we 
have used the tetraquark currents of type M to perform QCD sum rule analyses and investigated the newly observed 
ai(1420) state. 






















FIG. 4: The mass calculated using the currents (upper figures) and = 'tp^{qqqq) (lower figures), is shown 

with respect to the threshold value so (left figures) for M% = 1.2 (dotted), 1.5 (solid) and 1.8 GeV^ (dashed), and with respect 
to the Borel mass Mb (right figures) for so = 2.8 (dotted), 3.0 (solid), and 3.2 GeV^ (dashed). The working region is 1.0 
GeV^ < M% < 1.8 GeV^. However, The mass cnrves decrease quickly with M% from 1.0 GeV^ to 1.2 GeV^. Therefore, we 
choose the new interval 1.2 GeV^ < M% < 1.8 GeV^ as our working region. 


Combining the results of Ref. [2l| and the results obtained in this paper, we found that a mass prediction around 
1.5-1.6 GeV is often obtained (with respect to the threshold value sq). This may be a reasonable result: there are 
two tti states, ai(1420) and ai(1640), whose masses are close to each other, so that if a current couples to both 
of them and one still uses a one-pole parametrization, a mass value between these two masses would be obtained 
for this single pole. However, in the absence of definitive phenomenological models it is impractical to develop a 
two-pole parametrization that differentiates ai(1420) from ai(1640). In this paper we have used another approach, 
i.e., finding a current mainly coupling to ai(1420). We have used tetraquark currents of type M to perform QCD 
sum rule analyses, and found that the current leads to a mass of 1.44 ± 0.08 GeV. The good agreement of this 
result with the experimental value suggests that this current couples to the ai(1420) state supporting a tetraquark 
interpretation. 

We note that the quark content usus of the current means that it does not have a definite value of isospin(i.e, it 
is neither isospin one nor isospin zero). The isovector tetraquark current and its isoscalar partner can be constructed 
by changing the quark contents to be {usus — dsds) and {usus + dsds), respectively. However, the same sum rule and 
mass prediction would be obtained for all these three currents under SU{2) isospin symmetry. As noted in Ref. 
there is another isoscalar state, /i(1420), which has been well established in experiments [1^. It strongly couples to 
KK*, and is likely to be the isoscalar partner of ai(1420). If this is the case, our analyses would support tetraquark 
explanations for both of them. 

To conclude this paper, we study the possible decay channels of oi (1420). To do this, we use the Firez transformation 
and change the current with quark contents {usus — dsds), i.e., '4’M,itJ.{usus — dsds): 

i^M.iti{usus - dsds) = rtfC75S&(ua7/xCs^ - Ubjf_,Csl) + u^Cjf_,Sb{uaj5‘Csb - Ub'ybCsl) 


(30) 
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into a combination of {qq){ss) and {qs){sq) currents: 

^ {UaUa){sb^fij5Sb) “ {Ual^.'ybUa){sbSb) - {dada){sbl^.IbSb) + {dal^.l5da)(sbSb) i 
^ {UaltiSa)iSbl5Ub) “ {u al5 S a) {hi )iUb) - {daltiSa){hl5db) + {dal5Sa){sblt^db) , (31) 

^ {Ual''lbSa){h<^^ivUb) “ {Ua<^^iuSa){hl''Ib^b) “ {daYl5Sa){h^),,ydb) + {dad^^Sa){hl''Ibdb) , 

^(99)ss ^ {Ual''Ua){Sbd^vlbSb) “ (UaO'^i/75Ma) (s67''sb) - {dal’' da){hdfj.ul5Sb) + {dad ^vlbda){sbl''Sb) , 
through 

- dsJs) = • (32) 

We note that and both contain one qq meson and one ss meson, while and both contain 

one qs meson and one sq meson, where q represents an up or down quark, and s represents a strange quark. This 
equation suggests that ai(1420) may naively fall apart to: 

: ai(1420) ^ 0+ (ct( 600), ao(980), /o( 980) • • •) + 1+ (&i(1235), ai(1260) • • •) , 

: ai (1420) ^1-(X* (892) ■••)+0" (JiT •••) , (33) 

: oi(1420) ^ l+(i4:i(1270) •••) +1+(-fi^i(1270) •••) , 

4”^"" : ai(1420)^l-(p(770),w(782),(;i(1020)---) + l”(p(770),a;(782),<))(1020)---) . 

These are all iS-wave decay channels, while the possible P-wave decay channels can be obtained by naively relating 
dltilsq and 

4”^" : ai(1420) ^ 0+ (ct( 600), ao(980), /o( 980) • • ■) + 0" (^, r/, r?' • ■ •) , (34) 

4'^''^®'* : ai(1420) ^0-(P:---) + 1+(7^1(1270) •••) , (35) 

One extra constraint is that the final states of ai(1420) should contain one ss pair. Then the kinematically allowed 
decay channels are ^-wave ai(1420) —>■ K*{892)K and P-wave ai(1420) —>■ cr(600)?7, ai(1420) —> ao(980)7r and 
ai(1420) —> /o(980)7r. However, the P-wave decay channel ai(1420) —>■ d{600)g is forbidden by the conservation of 
isospin symmetry, and the P-wave decay channel ai(1420) —>• ao(980)7r is forbidden by carefully checking the detailed 
expression of . 

Summarizing all the above constrains, the possible decay patterns of ai(1420) are S'-wave ai(1420) —>• K*{892)K 
and P-wave ai(1420) —>■ /o(980)7r, the latter of which is observed by the COMPASS experiment Q. Similarly, 

the possible decay patterns of /i(1420) can also be studied, and they are S'-wave /i(142^ —AT* (892)AT and P- 

wave /i(1420) —)■ ao(980)7r, also consistent with the experiments Issl - f^ (see Refs. [38l - l40l| for related theoretical 
studies). 

BESIII is a good platform to carry out the search for ai(1420) by the J/ip radiative decay J/ip —>■ 7/o(980)7r 
if we take into consideration the strong interaction between ai(1420) and /o(980)7r indicated by the COMPASS 
measurement EMI- Of course, the production ratio of J/ip 7ai(1420) —> 7/o(980)7r is related to the inner structure 
of aiO420), which determines the strength of J/ip decaying into 7ai(1420). We note a former BESIII result in 
Ref. [3, where the J/ip —>■ 77r+7r“7r° and J/ip —>■ 77r°7r°7r° decays were studied and a large isospin violating process 
7^(1405) —>■ /o(980)7r° was observed [s^. If there exists a new state ai(1420), it may be revealed in a BESIII re-analysis 
of the 7^(1405) —i> /o(980)7r° branching ratio. This will provide a definitive test of whether 7^(1405) —>■ /o(980)7r° still 
has a large branching ratio when including the ai(I420) contribution in the J/ip —)• 77r+7r“7r° and J/ip —>■ 77r°7r°7r° 
decays. 


Appendix A: Other sum rules 


The sum rules using the currents = ip^^{qqqq) is 

= nM.i(so,M|) 


r( 


1 I (q^sGG) 

36864776 184327r6" 


'{gsqdGq)'^ (5sGG)(gg)" 


967r^ 


8647r2 


5 ( 99)4 , {qq){9sqdGq)\ _^^Mi^^ 

3671^ ^ 8772 4 

1 {9sGG){qq){gsqdGq) 

' Ml 576772 


(Al) 
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The sum rules using the currents = ip^^{qqqq) is 

= nM.2(so,M|) 

-n 


1 


, 184327r6 


{9sGG) 2 5{qq) ^ -s/mi 


^4 _ „2 


184327r6'' 


-( 


(gsqaGq)'^ {glGG){qqf 


487r2 8647r2 

The sum rules using the currents = ilj^^{qqqq) is 

= nM.3(so,M|) 

1 


+ 


187r2 47r2 

1 {glGG){qq){g,q(jGq) 


^ds 


Ml 


5767r^ 


/( 


122887r6' 


4 , {giGG) 2 5 {gg)^ {qq){gsq(jGq)\ _,/mI 

_1_ -^ — - c — - 1 ^ / S 


184327r6 


367r2 


s — 


87r2 






(gsqcrGq)'^ {glGG){qqY 


967r2 


2887r2 


1 {glGG){qq){gsqaGq) 
Ml 1927r2 


(A2) 


(A3) 


and the results are shown in Fig. [5] 



FIG. 5: The mass calculated using the current is shown with respect to the threshold value so (left panel) 

for Ml = 2.5 (dotted), 3.0 (solid) and 3.5 GeV^ (dashed), and with respect to the Borel mass Mb (right panel) for so = 2.8 
(dotted), 3.0 (solid), and 3.2 GeV^ (dashed). 


The sum rules using the currents = 'ilj^^{qqqq) is 


= nM,4(so,M|) 

r° ( 1 „4 I 5(gg)\ {m){9sq(TGq) \ _^/j^2^ 

Jq V 614471® 1843277® 187r2 47r2 / 

/ _ {gsqcrGq)^ _ {glGG) {qq)"^ \ 1 {g^GG){qq){gsqaGq) 

V 487r2 2887r2 ' Ml 1927r2 


(A4) 


and the results are shown in Fig. [6] 
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FIG. 6: The mass calculated using the current = 'il’^{qqqq) is shown with respect to the threshold value so (left panel) 
for M% = 2.5 (dotted), 3.0 (solid) and 3.5 GeV^ (dashed), and with respect to the Borel mass Mb (right panel) for so = 2.8 
(dotted), 3.0 (solid), and 3.2 GeV^ (dashed). 


The sum rules using the currents r]^ = ipi^^{qqqq) is 

= nM.5(so,M|) 


(A5) 


fSo 

J 4 m^ 


1 


368647r6 


s^- 


9607r6 


( {glGG) 7ms{qq) , ms(ss )^^2 


V 184327r6 
.2/^2 


3847r4 


1287r4 


^ 5{qq){ss) _ BmsjgsqcrGq) _ 13mf (g^GG) ^ ^ ^ {qq){gsSaGs) (ss){gsqaGq) 


V 3671^ 

/„2 


1927r4 


368647r6 

.2/;T„\2 


IGtt^ 


+ 


5127r4 15367r4 87r2 

■{gsqaGq){g,saGs) {glGG){qq)'^ {glGG){qq){aa) {glGG){33)‘^ 


167r2 

m^(ss)^ 


487r2 


Ams{qq)'^{33) 


V 967r2 34567r2 5767r2 34567r2 9 

ms{qq){33)'^ ms{glGG){gsqaGq) ms{glGG){gs3aG3) ml{qq){gsq(jGq) 


9 30727r4 

m'i{33){gsqcrGq) ml{qq){gs3aG3) 


167r2 


247r2 


)+m( 


92167r4 127r2 

1 ^ {glGG){qq){gsqaGq) 


23047r2 


{9‘iGG){qq){gs3oG3) {g‘^GG){33){gsqcrGq) {g'^GG){33){gs3aG3) 


7687r2 


7687r2 


23047r2 


msiqqrigsaaGa) 2ms(gg)(ss)(gs qcrGq) _ ms(gg)(ss}(gsSO'Gs) 

9 9 “ 12 

ms{33)'^{gsqaGq) mUgsqcrGq)'^ ml{gsqcrGq){gs3aG3) 


12 


967r2 


327r2 


)■ 
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The sum rules using the currents = fp^^{qsqs) is 

= UmAso,MI) 


4m2 


1 „4 „3 1 7 

(g^sGG) 

7ms{qq) mg{ss)\ 

1843277®" 48077®" ' V 

1843277® 

192774 64774 ) 


/ 5{qq){ss) _ 5ms{gsq(jGq) _ 17m;{ggGG) \ 

V IStt^ 9071^ 3686471® / 

{qq){gsSaGs) {ss){gsqaGq) ms{glGG){qq) bms{glGG){ss) ml{qqf 


Stt^ 


Stt^ 


‘iml{qq){ss) ^ ml{ss)'^ 


5{g‘iGG){qqf 


247r2 




5127r4 1530774 

{gsqaGq){gsSaGs) 


Stt^ 


487r2 


{g‘^GG){qq){ss) 5{g'^GG){ss)'^ ms{glGG){gsqaGq) 


34567r2 5767r2 34567r2 

5ms{ggGG){gsSaGs) 8ms{qq)‘^ (ss) 2ms{qq){ss)'^ 
92167r4 9 9 


3072774 

mg{qq){gsq<jGq) 

6772 


'ml{qq){gsS(TGs) _ ml{ss){gsqaGq) \ 1 / 5{g^GG){qq){gsqaGq) 
12772 g ^2 ) ^ 2304772 

{ 9 lGG){qq){gsSaGs) {g‘^GG){ss){gsqaGq) 5{g‘^GG){ss){gsSaGs) 


768772 


768772 


2304772 


2ms{qq)^{gsSaGs) Ams{qq){ss){gsqaGq) ms{qq){ss){gsSaGs) 


+ 


6 


ms{ss)‘^{gsqaGq) mUggqaGq)'^ ml{gsqaGq){gsSaGs) \ 


6 


48772 


16772 


)■ 


(A6) 
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The sum rules using the currents = 'ip^^^qsqs) is 

= nM.7(so,M|) 


rsQ 

J 4m2 


1 




11m, 


122887r6 25607r6 

5(gg)^ , 5(gg)(ss) 5(ss)2 


f {giGG) 7ms{qq) 23ms{ss) 


V18432776 3847r4 38477^ J 

5ms{gsgo'Gq) 5ms{gsScrGs) 23nn1{glGG) 


s — 


SGtt^ SGtt^ SGtt^ 
{qq){gsqaGq) {qq){gsSaGs) 


+ 


1927r4 

{ss){gsqcrGq) 


90774 


877^ 1077^ 1077^ 

ss){gsSaGs) 3ms{glGG){qq) ms{glGG){ss) ml{qq)'^ 


877^ 

Tnl{ssf 

1077^ 


512774 


512774 


3080477® 


^ml{qq){ss) 

877^ 


ds+ 


{gsqaGq}^ {glGG){qqf {g'^GG){qq){ss) 


9077^ 


1152772 


192772 


{9sGG){ss)^ {gsqcTGq){gsSaGs) {g^saGs)'^ ms{glGG){gsq(jGq) 


1152772 90772 90772 

ms{glGG){gssaGs) _ 14 ms(gg)^{ss) ms{qq){ss)‘^ 

^ 3072774 9 9 

ml{qq){gsSaGs) _ m'^^{ss){gsqaGq) \ 

24772 10772 


1024774 

5ml{qq){gsqaGq) 

12772 




+ 


1 ({glGG){qq){gsq(jGq) {g'^GG){qq){gsSaGs) {g'^GG){ss){gsq<jGq) 


Ml 


708772 250772 

{glGG){ss){gsS<7Gs) ms{qqf{gsSaGs) 

708772 ”^3 

ms{qq){ss){gsSaGs) ms{ss)‘^ {gsqaGq) ml{g‘^GG){ss)^ 


250772 

ms{qq){ss){gsqaGq) 

^m’j.igsqcrGq)'^ 


12 

■ml{gsqaGq){gsSaGs) 

32772 


12 


1152772 


32772 


and the results are shown in Fig. [T] 



(A7) 


FIG. 7: The mass calculated using the current = 'tp^^{qsqs) is shown with respect to the threshold value so (left panel) 
for Ml = 2.5 (dotted), 3.0 (solid) and 3.5 GeV^ (dashed), and with respect to the Borel mass Mb (right panel) for so = 2.8 
(dotted), 3.0 (solid), and 3.2 GeV^ (dashed). 
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The sum rules using the currents = ip^^{qsqs) is 

= UmAso,MI) 


(A8) 



1 „4 11^? „3 1 

(llig^sGG) 

7ms{qq) 

23ms{ss)\ 


614477® 128071® 

1 1843271® 

1927r4 

^ 1927r4 / 


+ - 


5{qq)^ 5{qq){ss) 5{ssy 5ms{gsqo'Gq) 5ms{gsScrGs) 163mg{ggGG) 


IStt^ IStt^ IStt^ 9671^ 

{m){gsqcrGq) {qq){gsSaGs) {ss){gsq(7Gq) 


487r4 


368647r6 


47r^ 


+ 


Stt'^ 


+ 


Stt'^ 


+ 


+ 




{ss){gsSa Gs) 3ms{g^GG){qq) ^ 5ms{gsGG){ss) ^ 2ml{qq)^ 3ml{qq){ss) 
47r^ 
m^(ss)^ 

87r2 

5{g^,GG){ss)^ 


{gsqaGqf _ 5{g^GG){qq)^ {g^,GG){qq){ss) 

487r2 11527r2 1927r2 

{gsqaGq){gsSaGs) (gsSaGs)'^ ms{glGG){gsq(jGq) 


11527r2 487r2 487r2 10247r4 

bms{glGG){gsS(jGs) 28ms{qq)'^{ss) 2ms{qq){ss)‘^ bml{qq){gsq(jGq) 


30727r4 9 

'ml{m){gsS(TGs) ■ml{ss){gsqaGq)' 


9 




127r2 


87r2 


+ 


+ 


1 f^{ 9 sGG){qq){gsq(TGq) {g'^GG){qq){gsSaGs) {g'^GG){ss){gsqaGq) 


Ml 


7687r2 2567r2 

5{g'^GG){ss){gsSaGs) 2ms{qq)'^{gsSaGs) 

7687r2 3 

mi,{qq){ss){gsSaGs) ms{ss)‘^ {gsqcrGq) 5ml{g'^GG){ss)‘^ 


2567r2 

2ms{qq){ss){gsqaGq) 

3ml{gsqaGq)'^ 


6 

mligsqcrGq) {gsSaGs) 


6 


11527r2 


167r2 


167r2 


and the results are shown in Fig. [S] 



FIG. 8: The mass calculated using the current = tjj^lqsqs), with respect to the threshold value so (left) for Ml = 2.5 
(dotted), 3.0 (solid) and 3.5 GeV^ (dashed), and with respect to the Borel mass Mb (right) for for so = 2.8 (dotted), 3.0 (solid), 
and 3.2 GeV^ (dashed). 
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